Logarithmic roughening in a growth process with edge evaporation 
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Roughening transitions are often characterized by unusual scaling properties. As an example 
we investigate the roughening transition in a solid-on-solid growth process with edge evaporation 
[Phys. Rev. Lett. 76, 2746 (1996)], where the interface is known to roughen logarithmically with 
time. Performing high- precision simulations we find appropriate scaling forms for various quantities. 
Moreover we present a simple approximation explaining why the interface roughens logarithmically. 



I. INTRODUCTION 



l-q 



A large variety of models for surface growth is known to 
display universal scaling laws [|l], ||, ^, Q ^ . In most cases 
one observes simple power- law scaling, i.e., the width W 
of an initially flat interface grows with time as 



(1) 



where 7 is the so-called growth exponent. In a finite 
system the width eventually saturates at a finite value 
Wsat, which grows with the linear system size L as 



(2) 



where a is the roughness exponent. Both asymptotic 
power laws can be combined in a single scaling form 



WiL,t)^L''f{t/L^), 



(3) 



where z = a/7 is the dynamic exponent and / is a scaling 
function with an appropriate asymptotic behavior. This 
type of power-law scaling, also known as Family- Vicsek 
scaling |6|, is generic for a large number of self-similar 
growth processes. However, in some cases the interface 
was found to roughen logarithmically with time. This 
may happen, for example, when a system undergoes a 
roughening transition. The purpose of this paper is to 
study a simple model with such a logarithmic scaling be- 
havior in more detail. To this end we consider a class 
of solid-on-solid growth models introduced a few years 
ago by Alon et al. |8| which have the special feature 
that atoms can only evaporate at the edges of terraces. 
The models exhibit a robust phase transition from a non- 
moving smooth phase to a moving phase, where the in- 
terface roughens continuously. Remarkably, this rough- 
ening transition takes place even in one spatial dimen- 
sion. At the critical point the interface width was found 
to increase logarithmically with time. However, so far 
it was impossible to unfold the complete scaling picture 
since previous numerical simulations were not accurate 
enough [^. Here we close this gap by performing high- 
precision simulations based on a parallelized code. As a 
main result wc find that Eq. (^) has to be replaced by a 
logarithmic scaling law of the form 



W'^{L,t) ^ XlnL + f{t/L'') 



(4) 



Moreover, we present a simple approximation which ex- 
plains why the interface roughens logarithmically. 
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FIG. 1: Dynamic rules for deposition and evaporation 
(^) of the restricted model introduced in Ref. Q. Note that 
evaporation from the middle of plateaus is not allowed. 



II. GROWTH MODEL WITH EVAPORATION 
AT THE EDGES OF PLATEAUS 



The growth models investigated in Ref. [[7| are defined 
as solid-on-solid deposition-evaporation processes with 
the special property that atoms can only evaporate at 
the edges of terraces. The models are defined on a d- 
dimensional lattice with L sites and periodic boundary 
conditions, where each site i is associated with an inte- 
ger variable hi. In Ref. ^ two variants of models have 
been considered. Here we will focus on the physically 
motivated restricted model in one dimension, where the 
heights at neighboring sites obey the inequality 



H+l 



< 1. 



(5) 



The model evolves by random-sequential updates, i.e., 
a site i is randomly selected and one of the following 
moves is carried out (see Fig. |^). With probability q an 
atom is deposited at site i, increasing hi hy 1. Otherwise 
one of the nearest neighbors is randomly selected. If the 
selected neighbor is at a lower height, indicating the edge 
of a terrace, an atom is removed from site i. In both 
cases a move is abandoned if the resulting configuration 
would violate the condition (^). Each attempted update 
corresponds to an average time increment At — l/L. 

In the limit L 00 the model has the following phe- 
nomcnological properties. If the growth rate q is small, 
the interface is flat and pinned to a spontaneously se- 
lected bottom layer. Small islands will occasionally grow 



2 



1.15 
1.1 
1.05 

'p 1 
0.95 
0.9 



0.85 









- 








q— U.ioojy 






q=0. 18860 






q=0. 18861 






q=0. 18862 






q=0. 18863 






q=0. 18864 











10 



10 



10 



10" 10 
t [MCS] 



10 



10 



FIG. 2: Density of exposed sites at the bottom layer no{t) 
times in a system with 32768 sites averaged over up to 1000 
independent runs for different values of q. 



on top of this bottom layer but will quickly be eliminated 
by desorption at the islands edges so that the system 
eventually reaches a stationary state characterized by a 
finite width. As q is increased, more and more islands 
on top of the bottom layer are produced, until above qc, 
the critical value of q, the islands merge and the interface 
starts moving. Approaching qc from below, the station- 
ary width is found to diverge logarithmically. 

In Refs. 0, 1 it was shown that the roughening transi- 
tion in this class of models is driven by a directed perco- 
lation (DP) transition Q at the spontaneously selected 
bottom layer. This means that the density of exposed 
sites at the bottom layer no can be interpreted as the 
density of active sites in a DP process. Therefore, at 
criticality this density is expected to decay as 



HQ '^t 



(6) 



where 6 = 13/ v± ~ 0.1595 is one of the dynamic critical 
exponents of DP in l-fl dimensions. Similarly the criti- 
cal behavior at the first few layers above the bottom layer 
can be described in terms of unidirectionally coupled DP 
processes Q . Roughly speaking each layer can be asso- 
ciated with a DP process which is coupled to the layer 
below. Effective couplings in opposite direction, which 
certainly exist in the model defined above, turn out to be 
irrelevant. Field-theoretic renormalization group studies 
revealed that the scaling exponents vw and (and hence 
their ratio z = v\\lv^^) are the same at all levels, while the 
exponents associated with the order parameters decrease 
with increasing height above the bottom layer. Thus the 
concept of unidirectionally coupled DP successfully ex- 
plains the critical behavior at the first few layers above 
the bottom layer. However, it cannot predict the criti- 
cal properties of the interface as a whole, especially the 
scaling of the interface width. 
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FIG. 3: The squared interface width and the average 
height H as functions of time in a log-lin plot. The simulation 
parameters are the same as in Fig. H. The inset shows that 
both quantities increase logarithmically at the critical point. 
Small deviations for t > lO'' can be traced back to the onset 
of finite size effects. 



III. NUMERICAL SIMULATIONS 



A. Scaling properties at the critical point 



Performing high-precision simulations on a parallel 

(I) to estimate the critical 
2 we obtain the result 



computer we first used Eq. 
threshold. As shown in Fig. 



9c 



1.8861(1) 



(7) 



which is much more accurate than the previous estimate 
reported in Ref. Q. Moreover, our simulations confirm 
that the dynamics at the bottom layer is indeed driven 
by a DP process. 



At the critical point q = qc the average height 



and the squared interface width 



1 ^ 

i=l 



(h, - H) 



(8) 



(9) 



are found to increase logarithmically with time as 

W^{t) ~ Tlni, (10) 



H{t) ~ Clint, 



(11) 



where r = 0.102(3) and a = 0.133(3) (see Fig. |). Thus, 
apart from different prefactors, the squared width and 
the mean height show the same type of logarithmic be- 
havior. We note that in a previous work the width was er- 
roneously conjectured to increase as W{t) ^ (Int)" "*^ [0. 
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FIG. 4: Finite size simulations: Squared width W^{L,t) and 
mean height H{L,t) averaged over at least 500 independent 
runs in a finite system with L — 8, 16, 32, . . . , 1024 sites for 
(a) periodic and (b) fixed boundary conditions. In the lower 
panel the curves for H{t) have been shifted vertically. 



B. Finite-size simulations 



In a finite system the critical behavior depends on 
the type of the boundary condition. In what follows we 
consider (a) periodic and (b) Dirichlet boundary condi- 
tions, where the sites at the boundary are fixed at zero 
height. As shown in Fig. ^, the average height H{L,t) 
and the squared width W^{L, t) first increase logarithmi- 
cally until finite-size effects become relevant and the sys- 
tem crosses over to a different regime, where the width 
saturates. Measuring the saturation levels we find that 



Apj -f \pj InL , 



(12) 



where the indices p, f stand for periodic and fixed bound- 
ary conditions, respectively. Our best estimates are: 



An 



0.161(2) 
0.163(2) 



(13) 



A; = -0.21(1), A/ 

Apparently the amplitude A is universal in the sense that 
it does not depend on the choice of the boundary condi- 
tion. 

Unlike the width, the interface height does only satu- 
rate if fixed boundary conditions are used. Here we find 
a similar formula 



where 



i7sat(i) ^S/+/i/lnL, (14) 



B/ = -0.24(1), /i/ = 0.212(10) . (15) 




FIG. 5: Finite size simulations: Data collapse of the curves 
shown in Fig. ^ according to the scaling forms ( [l7| ) for L = 
32, 64, . . . , 1024. In order to eliminate deviations due to initial 
transients all data points for t < 100 have been discarded. 



In order to collapse the curves shown in Fig. |j, recall 
that the roughening transition in this model is driven 
by a l-|-l-dimensional DP process at the bottom layer, 
which is characterized by a dynamic exponent 



z = l/u/ V 



1.5807. 



(16) 



Moreover it has been shown that the critical dynamics at 
the first few layers is the same as in a sequence of uni- 
directionally coupled DP processes, which are all char- 
acterized by the same dynamic exponent z. Therefore, 
it is reasonable to assume that the roughening transition 
as a whole is characterized by the dynamic exponent of 
DP so that finite-size scaling functions can only depend 
on the scaling-invariant ratio t/L'^. Thus, the expected 
scaling forms read 



W^{L,t) ~ \\nL + f{t/L'), 
H{L,t) ~ ^lnL + g(i/i^), 



(17) 



where /(^) and g{S^) are scaling functions with the asymp- 
totic behavior 



no 



A 



■Ine, 



InC 



(18) 



for ^ and f{Ci A and g{Ci = B for ^ ^ oo. In 
order to verify these scaling forms we plot W^^(L,<) — 
A InL and H{L,t) — /ilnL versus t/L^. As shown in 
Fig. H, we obtain excellent data collapses. 



Off-critical simulations 



A similar situation is encountered in off-critical simu- 
lations. In the subcritical regime q < qc the width and 
the average height saturate at 



ffsat(e) 



D 



7/ In |e| 
- (hlnlel 



(19) 



where e = q ~ qc denotes the distance from criticality. 
Following the same arguments as in the previous subsec- 
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FIG. 6: Off-critical simulations for e = ±0.0001, ±0.0002, 
. . . , ±0.0128 on a chain with 32768 sites averaged over 200 
independent runs. The figure shows data collapses according 
to the scaling forms (^9|). In order to eliminate deviations 
due to initial transients all data points for t < 100 have been 
removed. 



tion, the conjectured scaling forms read 

W\e,t) = -7?ln|e|±F(t|erii), 
Hie,t) = -c^\n\e\ + G{t\ep), 



(20) 



where ly^^ ~ 1.7338 is the temporal scaling exponent of 
DP in 1+1 dimensions. These scaling forms are expected 
to hold not only below but also above critical! ty. In fact, 
plotting W'^{e,t) + 7j\n\e\ and H{e,t) + (j)ln\ e | versus 1 1 e | ''ii 
(see Fig. 0), the best data collapses are obtained for 



77 = 0.173(10), = 0.230(10), 
while the constants are estimated by 

C = -0.40(1), £1 = -0.54(1) . 

D. Relations among the amplitudes 



(21) 



(22) 



Obviously the scaling forms (|17|) and (^^ suggest that 
the exponentiated quantities exp(VF^) and exp{H) dis- 
play ordinary power-law scaling. This would imply that 
the results for finite-size and off-critical scaling can be 
combined in a single scaling form, e.g.. 



exp[M^2(e,L,t)] = FieL^/"^ ,t/ L"") , 
exp[i7(e,L,t)] = GieL^/"^ ,t/ L') . 



(23) 



This in turn would imply that the six amplitudes 
A, /X, T, cr, 77, cj) introduced above, which now play the role 
of critical exponents, are no longer independent. Only 
two of them, say 77 and (f), are independent while the oth- 
ers can be expressed as 



jj. = (j)/v±_ , 
r = 77/1^11 , 



(24) 



where v± ~ 1.0968 and v\\ ~ 1.7338. As shown in Table| 
these relations are compatible with the numerical results. 



quantity 


amplitude 


estimate 




T 


0.102(3) 


m 


a 


0.133(3) 


W^{L,t) 


Xp 


0.161(2) 


W\L,t) 


A/ 


0.163(2) 


H{L,t) 




0.208(10) 


H{L,t) 


M/ 


0.212(10) 


W'{e,t) 


V 


0.173(10) 


Hie,t) 


<P 


0.230(10) 



TABLE I: Summary of the estimated amplitudes. 
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FIG. 7: Approximation of the average height by considering 
islands of size £1,^2,^3, ■. . (see text). 



IV. SIMPLE APPROXIMATION 



Why does the interface roughen logarithmically at the 
transition? To answer this question we present a simple 
approximation for the average interface height H{L) in a 
finite l±l-dimensional system at criticality. We consider 
fixed boundary conditions, i.e., the interface is pinned 
to the bottom layer at the two boundary sites. In be- 
tween the interface may return to the bottom layer sev- 
eral times, dividing the chain into intervals of different 
lengths ^i,£2, • ■ as sketched in Fig. 0. The approxima- 
tion relies on the following assumptions: 

1. The probability distribution for the interval sizes i 
is the same as the size distribution of inactive is- 
lands in a DP process, which is known to decrease 

as P(£) -^^/''^^^ ^^-0.25^ 

2. As we are using fixed boundary conditions, we will 
assume that this probability distribution is cut off 
at the system size, i.e., P{t) = for ^ > L. 

3. The sites, where the interface touches the bottom 
layer, can be regarded as fixed boundary conditions 
for the dynamics taking place at the following layer. 

While the first two assumptions can be justified by nu- 
merical checks, the third assumption is crucial since the 
sites at the bottom layer evolve in time and can only 
be regarded as fixed in an approximate sense. However, 
using this approximation we can derive a recursion rela- 
tion since each island (for example the island encircled 
by the dotted line in Fig. 0) can be regarded as an in- 
dependent growth process taking place in a finite system 



5 



of size ^ with fixed boundary conditions at the following 
layer. Obviously, the average height of such an island 
will be 1 + H{£). In order to calculate the average height 
of the entire interface H{L), one has to add up all these 
contributions weighted by the island size £. Replacing 
sums by integrals, this leads to the integral equation 



H{L) 



lodiiPie) 



(25) 



Shifting the denominator to the l.h.s. and differentiating 
with respect to L we obtain 



dH{L) 



dL 



[ deep{i) = LP{L), 

Jo 



or equivalently 



In / d££P{£), 



(26) 



(27) 



dH{L) _ d 
dL ~ dL 

leading to the solution 



H{L) = const + In / d£lP{t} 

f (28) 
= const H In L . 

Thus the approximation correctly reproduces the loga- 
rithmic increase of the average height in a finite critical 
system. However, the predicted amplitude /i = l3/v±_ ~ 
0.25 differs from the measured value /it w 0.21. 

A similar argument can be used to explain the loga- 
rithmic increase of the squared width. To this end we 
first calculate the second moment of the height 



M2(L) 



(29) 



Using the same assumptions as before, each island of 
size £ contributes to M2{L) with a term 



M2{£) 



= l + 2H{£) + M2{£). 



(30) 



Replacing sums by integrals we obtain the integral equa- 
tion 



M2{L) 



d££P{£) [1 + 2H{£) + M2{£)\ 



Iod££P{i) 
leading to the differential equation 
dM2iL) 



(31) 



dL 



d££P{£) ^ LP{L) + 2LP{L)H{L) . (32) 



In order to remove the mixed term on the r.h.s. we dif- 
ferentiate the squared width W^{L) = A'hiL) ~ H^L) 
with respect to L and combine the result with Eq. (|2q): 



d[W^L) - M2{L)] 
dL 



d£iP{£) = -2LP{L)H{L). 

(33) 



Inserting Eq. ( p^ ) the mixed term drops out and we ob- 
tain the differential equation 



dW^{L) 
dL 



d££P{£) = LP{L), 



(34) 



which has exactly the same form as Eq. (26). Therefore, 
VF^(L) and H{L) will only differ by a constant. Again 
the approximation correctly reproduces the logarithmic 
increase but cannot predict the precise value of the am- 
plitude. 



V. CONCLUSIONS 

In the present paper we have analyzed the scaling be- 
havior at the roughening transition of a restricted solid- 
on-solid growth process with evaporation at the edges 
of terraces. Restricting to the case of a one-dimensional 
substrate we have carried out extensive Monte Carlo sim- 
ulations, exceeding the temporal range of previous nu- 
merical studies by three decades. The numerical results 
confirm that the transition is indeed driven by a directed 
percolation process at the bottom layer. Furthermore, 
they clearly show that squared interface width W'^{t) 
and the average height above the bottom layer H(t) in- 
crease logarithmically with time. Analyzing numerical 
data from finite-size and off-critical simulations we have 
postulated appropriate scaling forms which generate ex- 
cellent data collapses. These scaling forms can be inter- 
preted in such a way that the exponentiated quantities 
exp(W^^) and exp(7J) obey ordinary power-law scaling, 
which allows one to derive various relations among the 
amplitudes. Moreover, we have presented a simple ap- 
proximation which explains why and H grow loga- 
rithmically. 

In contrast to the expectations expressed in we 
cannot find evidence for several competing length scales. 
Moreover, there is no evidence for logarithmic scaling 
based on local scale invariance ||ll[, rather the system 
shows ordinary scaling after an appropriate redefinition 
(exponentiation) of the order parameters. 

Our results suggest that any roughening transition, 
which is driven by a phase transition from an active into 
an absorbing state at the bottom layer, should exhibit 
this type of logarithmic roughening. An important ex- 
ample is a recently introduced model for dimer adsorp- 
tion and desorption [Q, where the transition is driven 
by a parity-conserving transition at the bottom layer. 
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